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We point out that there is no cosmological gravitino problem in a certain class 
of gauge-mediated supersymmetry-breaking (GMSB) models. The constant term in 
the superpotential naturally causes small mixings between the standard-model and 
messenger fields, which give rise to late-time decays of the lightest messenger fields. 
This decay provides an exquisite amount of entropy, which dilutes the thermal relics 
of the gravitinos down to just the observed mass density of the dark matter. This 
remarkable phenomenon takes place naturally, irrespective of the gravitino mass and 
the reheating temperature of inflation, once the gravitinos and messenger fields are 
thermalized in the early Universe. In this class of GMSB models, there is no strict 
upper bound on the reheating temperature of inflation, which makes the standard 
thermal leptogenesis the most attractive candidate for the origin of the observed 
baryon asymmetry in the present Universe. 



On leave from University of Tokyo. 



1 Introduction 



The minimal supersymmetric standard model (MSSM) is the most promising candidate 
for physics beyond the Standard Model (SM), since it naturally solves the "hierarchy 
problem" and leads to a successful unification of the gauge coupling constants [1] . Because 
SUSY is not observed in the real world, it should be broken around the TeV scale. Once we 
allow generic soft SUSY-breaking terms, we must face hundreds of new parameters, which 
makes the rate of flavour- changing neutral-current (FCNC) interactions many orders of 
magnitude larger than the present experimental bounds. In order to obtain a successful 
low-energy effective theory, various mediation mechanisms of SUSY-breaking effects have 
been proposed. 

The scenario most commonly considered in phenomenology is the minimal gravity- 
mediated SUSY-breaking (mSUGRA) models. This scenario is very simple and aestheti- 
cally attractive, since gravity does exist in nature. In the mSUGRA models, we also have 
a promising candidate for dark matter. That is the lightest SUSY particle (LSP), which 
is usually the lightest neutralino. The severest difficulty in this scenario is the lack of 
a natural explanation for the suppression of the FCNC interactions. A specific form of 
the soft SUSY-breaking masses must be imposed by hand in order to suppress the FCNC 
interactions. 

It has been argued that the suppression of the FCNC interactions can be naturally 
obtained in brane-world SUSY-breaking scenarios, such as anomaly- [2] and gaugino- [3] 
mediated SUSY-breaking models. In these scenarios, the fields relevant to SUSY breaking 
are assumed to reside on the hidden brane, which is geometrically separated from the 
visible brane where the SM fields are localized. Recently, however, a crucial observation 
has been made in Ref. [4], where the authors have found that the separation of the visible 
and hidden branes in a higher-dimensional space-time is not sufficient for suppressing 
the FCNC interactions. Consequently, we need additional ad hoc assumptions in those 
models. 1 

At present, the most attractive scenario seems to be the gauge-mediated SUSY break- 
ing (GMSB) [6]. In the GMSB models, the suppression of the FCNC interactions is 
realized in an automatic way, just because SUSY breaking occurs at a very low-energy 
scale. Furthermore, a whole spectrum of the superparticles in the MSSM sector is com- 
1 There is an attempt to realize anomaly-mediation models in a four-dimensional framework [5]. 
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pletely determined by only a few parameters, which allows us to discriminate the GMSB 
models from other candidates in the future collider experiments. 

Unfortunately, from a cosmological perspective, there exists a big difficulty in the 
GMSB models. That is the so-called cosmological gravitino problem. In these models, 
we have no natural explanation for the dark matter in the present Universe. Thermal 
relics of the gravitino, the LSP in the GMSB models, overdose the Universe once they are 
thermalized in the early Universe. In order to avoid the overproduction of the gravitinos, 
there is a severe upper bound on the reheating temperature of inflation Tr, which is 
about Tr ^5 10 6 GeV for 7113/2 = 10 MeV for instance, and it even reaches Tr^$ 10 3 GeV 
in the case of the lighter gravitinos m 3 / 2 ^ 100 keV [7]. Furthermore, we have to fine- 
tune the reheating temperature just below this upper bound to explain the required mass 
density of the dark matter. We also have to generate the observed baryon asymmetry at 
just the same reheating temperature. Therefore, for successful cosmology in the GMSB 
models, we need incredible fine-tunings of various parameters, which apparently belong 
to independent physics, such as SUSY breaking, inflation and baryo/leptogenesis. 2 This 
is a big drawback of the GMSB models with respect to the standard mSUGRA scenario. 

In this letter, we point out that there are indeed no such difficulties in a certain class 
of GMSB models. We consider direct gauge-mediation models where the SUSY-breaking 
effects are directly transmitted to the messenger sector without loop suppressions [9]. In 
this class of models, the specific form of the superpotential of the messenger sector is 
usually provided by the i?-symmetry: since it is violated by the constant term in the 
superpotential (W), which is anyhow required to cancel the cosmological constant, it 
is quite natural to expect that there are small mixings between the SM and messenger 
multiplets induced by the condensation (W). As a result, the lightest messenger particle 
decays into the SM particle and gaugino through the mixings. As we will see, the resultant 
late-time decays of the lightest messengers provide an exquisite amount of entropy, which 
dilutes the thermal relics of the gravitinos down to just the observed mass density of the 
dark matter in the present Universe. 

Surprisingly, this miracle turns out to be true almost irrespective of the mass of the 
gravitino and the reheating temperature, once the gravitinos and messenger particles are 

2 If there exist extra matter multiplets of a SUSY-invariant mass of the order of the "/x-term", the 
observed baryon asymmetry and gravitino dark matter can be simultaneously explained in a way totally 
independent of the reheating temperature [8] . 
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thermalized in the early Universe. Consequently, by assuming natural mixings between 
the SM and messenger fields, the severe upper bound on the reheating temperature is 
completely eluded. This fact makes it much easier to construct a realistic inflationary 
scenario in the GMSB models. This result also has an important implication on the origin 
of the observed baryon asymmetry. We will see that the standard thermal leptogenesis, 
through out-of-equilibrium decays of the right-handed Majorana neutrinos [10], is now 
the most attractive mechanism to generate the observed baryon asymmetry in the GMSB 
models. 



2 Required Amount of Entropy 

In the GMSB models, the longitudinal component of the gravitino (~ Goldstino) ip in- 
teracts fairly strongly with the SM particles. The total production cross section of ip by 
scattering processes is given by [7] 

(S scatt t; rc i) « 5.9 ^ , (1) 

where ( ) denotes the thermal average; M* = 2.4 x 10 18 GeV is the reduced Planck scale, 
mg is the mass of the gluino and g 3 is the coupling constant of the SU(3)c gauge group. 
Then, the resultant interaction rate is given by 

1 scatt ~ (S sca ttf re l) ^rad j (2) 

with n rad = (C(3)/7r 2 ) T 3 being the number density for one massless degree of freedom. 

By the scattering interactions, the gravitinos are kept in thermal equilibrium if T SCSutt /H 
^ 1, where H is the Hubble parameter of the expanding Universe. The corresponding 
freeze-out temperature of ip is estimated to 

T , « i Tev (tm\ l,2 f^y njnvY _ (3) 

' V 230 J U0 keV/ V m £ J 

where g*(Tj) denotes the effective massless degrees of freedom when the cosmic tem- 
perature T = Tf. 3 If the reheating temperature of inflation T R is higher than Tf, 

3 For the light gravitino m 3 / 2 ^S 10 keV, the decay processes of SUSY particles are comparable with 
scattering ones and keep the gravitinos in thermal equilibrium until the temperature drops below the 
superparticle-mass scale. 
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gravitinos are in thermal equilibrium in the early Universe. Here, we have assumed the 
radiation-dominated Universe at the freeze-out time of the gravitinos. We will justify this 
assumption later, even in the presence of the messenger particles. 

If there is no additional entropy production, the resultant yield of the thermal gravitino 
is estimated to 45 3 3 

Y3/2 (= ^f) = 2^(T f ) — (2) ' (4) 

where 113/2 is the number density of gravitinos and s is the entropy density. In terms of 
the density parameter, it is written as 



Vl 3/2 h 2 ~ 5.0 x 



m^_\ ( J30\ 

lOkeVy U( T /)/ ' 



where h is the present Hubble parameter in units of 100 km sec -1 Mpc -1 , and fi 3 / 2 = 
P3/ 2 /p c ; P3/2 an d p c are the energy density of the gravitino and the critical density in 
the present Universe, respectively. Since the observed mass density of the dark matter 
is Qdm^ 2 — 0.1-0.2, the required dilution factor via the late-time entropy production is 
given by 

The above entropy should be supplied after the freeze-out time of the gravitinos. 



3 Decay of the Lightest Messenger and fi 3 / 2 

In this section, we discuss the decay of the lightest messenger and the resultant mass 
density of gravitino dark matter. Further constraints will be discussed in the next section. 
In this letter, we adopt the simplest messenger sector, which consists of a pair of chiral 
supermultiplets $ + <£>, describing a Dirac fermion of mass M and two complex scalar 
fields of mass squared M 2 ± F, where F is the F-term SUSY-breaking component of the 
mass of the messenger multiplets. In order to preserve the success of the gauge-coupling 
unification, we assume that $ + $ transform as 5 + 5 under the SU(5)gut gauge group. 4 

Under this setup, the gauginos and SUSY particles in the MSSM sector obtain the 
following soft SUSY-breaking masses, M a and m s 2 oft , via gauge interactions at the one- 
adopting a pair of 10 + 10 messenger multiplets does not change the basic results in this letter. 
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and two-loop level, respectively: 

M^A, roLt . 2? { Ca (_|L) 2 } A2 , ( 7 ) 

where g a (a — 1, 2,3) are gauge coupling constants and C a are the quadratic Casimir. 
Here, A = F/M determines the overall scale of the soft SUSY-breaking masses; A ~ 
10 5 GeV is required to obtain the correct size of soft masses. The mass of the messenger 
particles M and the gravitino mass m 3 / 2 are related by the formula: 



M 



v^Of* , in8 . / m 3/2 \ /10 5 GeV\ A7 



Here, we have defined k = F/Fdsb < 1) where Fdsb denotes the original F-term in the 
dynamical SUSY-breaking sector. In direct gauge-mediation models, this ratio is not 
loop-suppressed and naturally given by k^, 1. 

In the present work, we assume that there is the following mixing term between the 
SM and messenger multiplets through the small i?-symmetry-breaking effects caused by 
the constant term in the superpotential, (W): 5 

(W) 

SW = f \^5 M 5 = /m 3 / 2 5 M 5 , (9) 

where / denotes some unknown coefficient of the order of 1 and the subscript M denotes 
the multiplet that belongs to the messenger sector. In the present scenario, the lightest 
messenger superparticle is most likely the scalar component of a weak doublet. By virtue 
of the above small mixing term, the lightest messenger can decay into a SM lepton and a 
gaugino. 6 The decay rate is estimated to be 



16n \ M J 

The resultant decay temperature of the lightest messenger is given by 



(10) 




5 Hcre, we assume, for instance, the i?-charge for 5 and 5m to be +1 and — 1, respectively. 
6 We have neglected the decay channels through small Yukawa interactions. 
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Here, we have used the relation in Eq. (8). By comparing with Eq. (3), one can see 
that the decays of the lightest messenger always take place after the freeze-out of the 
gravitinos. Therefore, the entropy production associated with the decays of the lightest 
messengers dilutes the thermal relics of the gravitinos. 

Now, let us estimate the amount of entropy produced by the late-time decays of 
the lightest messengers. We will also justify the assumption of the radiation-dominated 
Universe made in deriving Eq. (3). Assuming the stability of the lightest messengers, 
their relic density has been estimated in Ref. [11] as 

M^H, 5X1 o-G^y, (») 

/ M \ 2 

^ 2 -10 5 (— — ) , (13) 



,10 6 GeV, 

where um/s denotes the frozen-out value of the yield of the lightest messengers and Qm 
is the corresponding density parameter. After the freeze-out of the lightest messengers, 
the total energy density of the Universe is given by 

7T 2 27T 2 

p = ^9*{T)T* + —g*(T)T 3 MY M , (14) 

where the first and second terms represent the contributions from the radiation and from 
the lightest messengers, respectively. The thermal relics of the lightest messengers begin 
to dominate the energy density at 

where we have used the relation in Eq. (8). From Eqs. (3) and (15), one can see that 
the matter-dominated Universe starts well after the freeze-out time of the gravitinos, and 
hence the assumption made in the derivation of Eq. (3) is justified. 

By assuming the instantaneous decays of the lightest messengers, which is accurate 
enough for the present purpose, we can obtain the dilution factor from energy conservation 
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as' 



( Saftcr \ _ 4 MYm 
V ^before / 3 Tfi 

~ 4.9 x 10 2 — — — , 16 

V10 8 GeV J V T d ) ' K ' 



7 Wc arc grateful to K.Hamaguchi for pointing out an error in the previous version. 
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where Sbefore (safter) denotes the entropy density of the Universe before (after) the decays 
of the lightest messengers. By substituting Eqs. (8) and (11) into Eq. (16), one can see 
that Am has the correct order of magnitude of the required dilution factor A in Eq. (6). 
From Eqs. (5), (8), (11) and (16), we can derive the final gravitino abundance as follows: 

1 



n 3/2 h 2 



n 3/2 h 2 



0.14 x / 




\ 5/2 

3 x 10 5 GeV J 



A/k 



2 keV 

, m 3 / 2 , 



1/2 



;i7) 



Astonishingly, a natural parameter k ~ 0.1 - 1 in direct gauge-mediation models and 
A fa 10 5 GeV, which is needed so as to obtain the correct size of the soft SUSY-breaking 
masses, leads to just the mass density of the gravitinos required to be the dominant 
component of dark matter in the present Universe. Furthermore, the resultant abundance 
of gravitinos has only a mild dependence on its mass mz/2- These facts can clearly be 
seen from Fig. 1, where we show a contour plot of fl 3 /2h 2 in a (7713/2 - k) plane. One 
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Figure 1: Contour plot of the abundance of the thermal gravitinos VL 3 /2h 2 in a (m 3 /2 - k) plane. 
The solid lines correspond to Q 3 /2h 2 — 0.3, 0.03 from the bottom up, respectively. Here, we 
have fixed other parameters as g*{Td) = 10, g*(Tf) = 230, A = 10 5 GeV, / = 1, for simplicity. 

can see that, for a given k-factor, the gravitino masses within one order of magnitude 
around a certain value give rise to a cosmologically interesting mass density of dark 
matter (0.03 ^ f2 3 / 2 ^ 2 ^ 0.3). As a result, the gravitino dark matter is a very natural 
consequence of the direct gauge-mediation models with the mixing term between the SM 
and messenger multiplets given in Eq. (9). 
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Before closing this section, we briefly comment on the property of the resultant grav- 
itino dark matter. The free-streaming length of the gravitino is given by [12] 

The late-time entropy production does not change the result in this approximation. For 
a more accurate estimation, see Ref. [13]. The gravitino of a mass in the range m 3 / 2 — 
(1-1.5) keV, whose free-streaming length is about Rf ~ 0.1 Mpc, is an interesting warm 
dark matter candidate, 8 which may reconcile the predictions of the cold dark matter with 
observations [15, 13]. Heavier gravitinos with mass m 3 / 2 > (a few) keV serve as the cold 
dark matter. The gravitino with mass m 3 /2 1 keV is now disfavoured by observations 
of Lyman-a forest [16] and the history of cosmological reionization [17]. 



4 Further Constraints 



In this section, we consider the non-thermal gravitino production from decays of the next- 
to-lightest superparticles (NLSPs). This contribution may spoil the successful prediction 
in the previous section, since a huge number of NLSPs are produced in the decays of the 
lightest messengers below their freeze-out temperature. In addition, if the decay process 
of the NLSP takes place when T ^ 5 MeV, it may also spoil the success of the Big Bang 
neucleosynthesis (BBN) [18]. 

Let us first consider the constraint from the BBN. The NLSP, which is the bino or 
stau in the GMSB models, decays into its superpartner and a gravitino, with the following 
decay width: 

1 m b x 

Tx ~ 48^ m\ l2 Ml ' (19) 
where m x is the mass of the NLSP. The corresponding decay temperature of the NLSP 
is given by 

T x = (™\ V4 Jfjf, * 5 MeV ( ^H 5 ' 2 (l™*) . (20) 

8 The possibility of the gravitino warm dark matter with a small entropy production has been discussed 
in Ref. [14], where the authors introduce the extra matter of mass Mx — 10 12 GeV to obtain the required 
lifetime of the lightest messenger particle. For their scenario to work, we need an ad hoc tuning on Mx- 
Therefore, their model does not solve the fine-tuning problem for the gravitino dark matter stressed in 
the introduction of this letter. 
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In order not to spoil the success of the BBN, the decays of the NLSPs should be completed 
before the onset of the BBN, T x ^ 5 MeV. If we take the natural range of the mass of 
the NLSP, for instance m x ^ 250 GeV, m 3 / 2 ^5 10 MeV is required for the success of the 
BBN. 

Let us now estimate the contribution of the non-thermal gravitinos to the total mass 
density of the dark matter. If Ym > F x , the produced NLSPs have enough time to 
annihilate before they decay into gravitinos. This case corresponds to 

•"^"^(lootv) ( A/ k ) ■ (21 > 

In this case, the resultant yield of the NLSPs before they decay is given by the following 
value, to a good approximation [19]: 



45 W " (22) 



87^(7^) M*T d ' 




where (crv) x denotes the s-wave annihilation cross section of the NLSP. The subsequent 
decays of the NLSPs produce the same number of gravitinos. Therefore, combined with 
Eq. (11), the resultant abundance of the non-thermal gravitinos is given by 

r-2\ 

• (23) 

In the case of the stau NLSP, the s-wave annihilation cross section is approximately 
given by (crv) x ~ 10~ 7 GcV~ 2 (100 GeV/m x ) 2 . Even in the case of the bino NLSP, 
(av ) x ^ 10~ n GeV -2 is naturally obtained for relatively large tan/3, which is about tan/5 
^ 15 (30) for m x ~ 100 (250) GeV, via small higgsino contamination. 

The above estimation is not valid if Y M < T x , since the produced NLSPs immediately 
decay into gravitinos. In this case, we have to follow the full evolution by solving coupled 
Boltzmann equations. 9 The result is given in Fig. 2. In this calculation, we have conser- 
vatively assumed that the number of NLSPs produced per decay of the lightest messenger 
particle is given by N x m M/E t h, where pa m x /Td is the threshold energy to produce 
the NLSPs by scatterings with the thermal backgrounds with T = T^. One can see that 
the contribution of the non-thermal gravitinos is smaller than several per cent of the total 
mass density of the dark matter as long as (av ) x ^ 10 -11 GeV -2 . 
Significant annihilations take place also in this case. 
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Figure 2: Contour plot of the abundance of the non-thermal gravitinos Q^^h 2 . Solid (dashed) 
lines correspond to f2?wf!/i 2 



&3/ 2 " — 10, 5, 2 and 1% of the total mass density of the dark matter 
( ^dm^ 2 — 0.15), for m x = 250 (100) GeV, from the bottom up, respectively. The previous 
estimation given in Eq. (23) is valid in the region where the contour lines are almost straight. 
In this calculation, we have fixed k = 1. 



In summary, the prediction of the gravitino dark matter given in Eq. (17) is not 
spoiled in most of the parameter region. The relevant constraint comes only from the 
BBN, which is satisfied as long as m 3 / 2 ^5 10 MeV for a reasonable range of m x . Finally, 
we briefly comment on another interesting aspect. The non-thermal gravitinos produced 
from the NLSPs behave as a hot/warm component of dark matter, which contributes 
several per cent to the total mass density of dark matter in the case of the bino NLSP 
with moderate tan/3. Therefore, the present model can naturally realize a mixed dark 
matter scenario. By virtue of the late-time production of the NLSPs, an unnaturally 
large hierarchy between the slepton and the bino masses is not required to realize a mixed 
dark matter scenario [20]. 



5 Thermal Leptogenesis 

As we have mentioned in the introduction, the present model has important implications 
on the origin of the baryon asymmetry in the present Universe. By virtue of the entropy 
production by the decays of the lightest messengers, we obtain the required mass density 
of the dark matter without any adjustments of the reheating temperature. Consequently, 
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the standard thermal leptogenesis [10, 21] now becomes the most promising candidate for 
the origin of the observed baryon asymmetry in the GMSB models. 

Let us start by introducing the relevant terms in the superpotential: 



W = -M Rl N t N t + h ia N t L a H u , 



(24) 



where JVj (i = 1,2,3) denote the heavy right-handed Majorana neutrinos of mass Mm', 
L a (a = e, /i, r) and H u denote the lepton doublets and the Higgs doublet that couples to 
up-type quarks, respectively The lepton-number asymmetry per decay of a right-handed 
neutrino JVj is given by [10, 22] 



E.r(jv, 



L a + H u ) - E a T(Ni ► L a + H u ) 



1 



L a + Hu) 



1 

'8tt 



^Elm[{(M%} 



v 



M 2 Ri . 



Ml, 
M 2 m , 



(25) 



where JVj, L a and H u (L a and H u ) symbolically denote fermionic or scalar components of 
corresponding supermultiplets (and their antiparticles); Tv{x) and J-s(x) represent the 
contributions from vertex and self-energy diagrams, respectively 



Fv{x) = In ^1 + , Fs(x) 



x — 1 



(26) 



For hierarchical right-handed neutrinos Mm <C M R2 ,M R3 , the lepton asymmetry is 
dominantly supplied by decays of the lightest right-handed (s)neutrinos. In the following 
discussion, we assume that this is the case for simplicity. In this case, the expression of 
the asymmetry parameter e\ is given by 



ei = 



3 M R1 mv3 



(27) 



where m„3 is the mass of the heaviest left-handed neutrino and 5 c g is an effective CP- 
violating phase. If the lightest right-handed (s)neutrinos are in thermal equilibrium, the 
resultant lepton asymmetry is given by the following formula: 



1 



x 



45 



C(3) 



A 2n*g*(T B ) vr 2 



(28) 
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where Tb is the freeze-out temperature of N\, and g^Ts) ~ 270; 10 k denotes the fraction 
of the produced asymmetry that survives washout processes by lepton-number-violating 
interactions after N\ decay. Here, we have assumed that there is no entropy production 
before the decays of the lightest messengers, and set Am = A for simplicity. For k ~ 1, 
the following out-of-equilibrium condition should be satisfied [21]: 

mi = ^F r7vi = {hh ' )n ^- 5 x 10-3 ev • (29) 

The lepton asymmetry produced by the right-handed (s)neutrino decays is subsequently 
converted into the baryon asymmetry by the sphaleron effects: 

^ = , (30) 

s s 

where C is a number of 0(1), which takes the value C = —8/23 in the MSSM. 

From Eqs. (6), (27), (28) and (30), we can derive the baryon asymmetry in the present 
Universe. In terms of the density parameter it is written as 

In Fig. 3, we show the lower bound of M R1 (and hence, it is the lower bound on the 
reheating temperature Tr) to obtain the required baryon asymmetry. Here, we have 
assumed the following relation for simplicity: 




k5 cS < 1. (32) 

As can be seen, the observed baryon asymmetry f2#/i 2 ~ 0.02 can be naturally explained 
by the thermal leptogenesis in a wide range of the gravitino mass. 11 From Eq. (8), one 
can see that the required lower mass bound on the right-handed neutrino is always larger 
than the mass of the messenger fields, and then the thermalization condition for the mes- 
senger particles is always satisfied. We should stress that we do not have to fine-tune the 



10 In this expression, we have assumed that a pair of messenger multiplcts (5 + 5) are in thermal 
equilibrium. After the decoupling of the messenger particles, the contents of the MSSM give rise to 
5* i=a 230. If we take into account the effective degrees of freedom in the dynamical SUSY-breaking 
sector, g*(Te) would be larger than this value, but it changes the resultant asymmetry by only a factor 
of at most 0(1). 

11 In the simplest chaotic inflationary scenario, for example, 0(1) couplings of the inflaton to the SM 
fields lead to T R w 10 13 GeV. 
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Figure 3: The lower bound on Mm (and hence Tr) to obtain the observed baryon asymmetry. 
Here, we have assumed the hierarchical right-handed neutrinos (Mm <C Mr2, Mrs) and the 
relation given in Eq. (32). The shaded region is disfavoured by the BBN constraint. 



couplings of the inflaton to the SM particles to reduce the reheating temperature, which 
is usually imposed in SUSY inflationary models because of the cosmological gravitino 
problem. We should also stress that there is an interesting possibility to determine 
directly in future collider experiments if the gravitino is lighter than about 100 keV [24]. 
Note that such a light gravitino has been considered as unlikely because of the cosmo- 
logical gravitino problem. However, as we have seen, it is now indeed favoured from 
cosmological perspectives. 



6 Conclusions and Discussion 



In this letter, we have pointed out that there is no fine-tuning problem to obtain the 
required mass density of the dark matter in a certain class of GMSB models. By virtue 
of the small mixing between the SM fields, which is induced by the i?-symmetry-breaking 
effects, the lightest messenger particle has a finite lifetime and provides an exquisite 
amount of entropy, which dilutes the thermal relics of the gravitinos down to just the mass 
density required for the dark matter. This phenomenon takes place naturally, regardless of 
the gravitino mass and the reheating temperature of inflation as long as the gravitinos and 
messenger fields are thermalized in the early Universe. There is no severe upper bound 
on the reheating temperature in this class of GMSB models, which makes the standard 
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thermal leptogenesis very attractive as the origin of the observed baryon asymmetry. 

The present scenario should have important implications also on other candidates for 
the origin of the present baryon asymmetry. It would be very interesting to reanalyse 
those models with the disappearance of the cosmological gravitino problem taken into 
account. 
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